Phase Diagram of a Lattice Microemulsion Model 
in Two Dimensions 

Per Arne Slotte 
Gruppe for teoretisk fysikk 
Institutt for fysikk — NTH 
N-7034 Trondheim 
Norway 

The phase diagram of a lattice microemulsion model proposed by Ciach, R~0ye and 
Stell [jlJQ is studied using mean- field theory and Monte Carlo simulations. Surfactant 
directional degrees of freedom are summed out exactly before mean-field theory is 
applied, and the resulting phase diagrams are much improved compared with previous 
mean-field results. The critical line and tricritical point is located using Monte Carlo 
simulations and finite size scaling. 
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I. INTRODUCTION 



Oil, water, and surfactant mixtures exhibit very special properties that makes them 
interesting from a practical as well as from a theoretical point of view. Under certain 
circumstances they form microemulsions — phases in which microscopic oil- and water- 
occupied regions are separated by thin layers of surfactant. The characteristic length scale 
of this structure is typically much larger than the "natural" length scale determined from 
the particle interactions, while still microscopic. 

This paper discusses the phase diagram of the two-dimensional version of a lattice model 
for a three-component microemulsion system proposed by Ciach, H_0ye and Stell and 
studied in subsequent papers H^HH^H • ^ ne f° now i n g I wm refer to this model as 
the CHS model. Figure |](a)-(c) shows the three types of phase diagram topologies that 
is expected for different surfactant strengths. In the following I will refer to these 

phase diagram topologies as type A, B and C respectively. The phase diagrams are given 
in the temperature/surfactant-chemical-potential plane for equal oil/water chemical poten- 
tials. With a weak surfactant, the phase diagram is of the simple three- component type 
shown in figure |](a), with two regions, disordered and oil/water coexistence, separated by 
a line of continuous transitions which changes into a first-order line at a tricritical point. 
When the surfactant is strong, there are two main possibilities shown in figure |I|(b) and (c). 
Both diagrams have four regions: disordered, microemulsion, oil/water coexistence, and or- 
dered structured phases (incommensurate, layered, bicontinuous . . . ). The disordered- and 
microemulsion regions are separated by a Lifshitz line PJlOU where the peak in the water- 
water structure function moves away from zero wave vector, and not by a phase transition. 
Note that some workers use the term disorder line for the Lifshitz line. I use here the nomen- 



clature of Gompper and Schick [|10j, which discriminates between the Lifshitz line, and the 
disorder line at which the asymptotic decay of the water-water correlation function changes 
from monotonic to nonmonotonic. In phase diagrams of type B, the disordered region is 
separated from the oil/water coexistence region by a line of continuous phase transitions. 
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The transition changes into first order at a tricritical point. The Lifshitz line intersects 
the first-order line, and thus we have the characteristic oil/water /microemulsion coexistence 
observed experimentally in many systems. The microemulsion is separated from the ordered 
structured phases by lines of phase transitions that may be either continuous or first order. 
In phase diagrams of type C, the line separating disordered from oil/ water coexistence is 
always a line of continuous phase transitions, and the line ends at a Lifshitz point where 
it intersects the Lifshitz line. Thus there is no oil/water /microemulsion coexistence. The 
microemulsion is separated from an ordered structured phase by a line of continuous transi- 
tions near the Lifshitz point, but the transition may be of any order elsewhere. The range of 
parameters where phase diagrams of type C is found is substantially reduced in the present 
calculations, compared to earlier mean- field calculations @J^|. 

The article is organized as follows: In section |H] the CHS hamiltonian is rewritten in 
terms of several spin-1 Ising variables on each site. In section |f| an exact mapping of 
the model onto an ordinary spin-1 model with temperature-dependent (effective) multispin 
couplings is developed. The connection between the CHS model and the spin-1 model of 
Schick and Shih |T3| is also discussed. In sections [TV] and [V] standard mean-field theory is 
applied to the effective spin-1 model. Results from Monte Carlo simulations are presented 



in section [VT|. In section [VII] the resulting phase diagrams are discussed and compared with 
phase diagrams obtained in related works. 



II. THE MODEL 

Each site in a hypercubic <i-dimensional lattice is occupied by oil, water or surfactant. 
The surfactant can take 2d different orientations along the lattice directions. This gives a 
total of 2 + 2d states for each site. 

The hamiltonian may be written in terms of 1 + d spin-1 Ising variables per site: a and 
r Xi , where x% is a lattice direction, o = ±1 represents oil/water, a = surfactant, and the 
r's different surfactant directions. Note that these variables are dependent since one, and 
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only one, of the variables at a given site will take a nonzero value. This gives again a total 
of 2 + 2d states. 

For simplicity it is assumed that the hamiltonian is symmetric with respect to interchang- 
ing oil with water and flipping the surfactant ends, ie a — ► —a and r — > —r. To be able 
to regulate the oil/water concentration ratio, the hamiltonian must in addition include an 
unsymmetric field term. If we include nearest-neighbor interactions only, and exclude direc- 
tion dependent surfactant-surfactant interactions, the most general form of the hamiltonian 
is: 

n = - £ ( H °i + ^) - E + K °l°) + A (W" - ' ) • (i) 

Here the first sum runs over all sites, while the second sum runs over all pairs (ij) of 
nearest-neighbor sites, and xjj is the lattice direction parallel with the bond (ij). The 
physical significance of the parameters in the hamiltonian are: H regulates the oil/ water 
concentration ratio. \x regulates the surfactant density. J > ensures that oil and water do 
not mix at low temperatures, and will be used as a temperature scale. K is the isotropic 
surfactant-surfactant interaction, I will assume K > 0. A is the surfactant strength, the 
sign of A is physically irrelevant. Note the antisymmetric form of the v4-term, this reflects 
the amphiphilic nature of the surfactant (ie. that one end loves water and the other end 
loves oil) . This term is the only term in the hamiltonian that is new compared with a simple 
hamiltonian for a three-component system with isotropic interactions. 

The hamiltonian has in other work [|T],^|,^|,^,^|J^J^] been formulated in terms of the oil-oil 
(water-water) energy, —b, the oil-surfactant energy, ±c, and the chemical potentials for 
oil, jLti, water, /i 2 , and surfactant, fi^ = \i± = ^ = fiQ = fj, s . The correspondence between 
the parameters in the hamiltonian (P above and the parameters b, c, fix, fi<i, and fig, is as 
follows: 

H = ^(/il ~ fl2) 

. 1. 



A = c (2) 
H = n s - -(/ii + // 2 ) 

if = -6 
2 

The ground states, and T = phase diagram, of the model have been previously analyzed 
for the K = J case ||, and for all K > — J the structure is the same. This is summarized 
in figure ||. 

The present spin-language formulation may not be appropriate if direction dependent 
surfactant-surfactant interactions are important, since the inclusion of such terms will com- 
plicate the hamiltonian ([]]) substantially 

It should be noted that the purpose of the present lattice model is not to accurately model 
the phase diagram of the real oil-water-surfactant system. A simple model of this type does 
not even model the two component (oil-water) system correctly: Due to the complexity of 
the oil-oil, water-water and oil-water interactions a temperature-dependent coupling J is 



needed to reproduce the important property called reentrant solubility, [|lT| , |i2|| i.e. that at 
fixed concentration oil an water mixes at high temperatures, and phase separates below a 
certain temperature, but mixes again in a range of lower temperatures. This stems from the 
complex nature of the bonds which can be in one of several van der Waals and hydrogen 
bonding states. Moreover, oil and water are certainly not symmetric in their surfactant 
interactions. The aim is rather to study the simplest possible model to be able to isolate the 
basic ingredients that are necessary to reproduce the key properties of real microemulsions. 

III. CORRESPONDENCE WITH SPIN-1 ISING MODEL 

The CHS model is closely related to two other lattice models proposed for microemulsion 
systems. One of these is the model proposed by Matsen and Sullivan 0,0] which is a 
generalization of the CHS model. The other is the spin-1 Ising model proposed by Schick 
and Shih [TB| and studied in a series of subsequent papers [|I^,|n],|T^.|T^.|T^.|E| . In the following 



I will call this model the Schick model. 



It is clear from the form of the hamiltonian (|l[) that the Schick model and the CHS model 
are closely related since both are extensions of the basic spin-1 Blume-Emery-Griffith (BEG) 
model: 



ftBEG = - E ( Hcr i + f^i) ~ E (JWi + K °l°1) ■ ( 3 ) 

» (*i> 

In the Schick model, the microemulsion character arises from an additional three particle 
interaction L<7j(l — er|)o"fc, where the sites ijk are neighbors on a line. The CHS model 
is more "realistic" with directed surfactants and pair interactions. As long as no direction 
dependent surfactant-surfactant interactions are included in the CHS model, the directional 
degrees of freedom can be summed out independently on each site resulting in an effective 
spin-1 hamiltonian of the Schick type with temperature-dependent multiple site couplings 



Ll|, |20|jn| . Hence, there exist an exact mapping between the CHS model and an extended 
version of the Schick model. Below, this mapping is explicitly calculated for the one dimen- 
sional version of the CHS model. The explicit expression for the effective hamiltonian in 
two dimensions is given in appendix A. 

The CHS hamiltonian ([!]) has the following symmetries that also must be valid for the 
effective spin-1 hamiltonian: 

• A global spin flip a — > —a, H — ► — H r — > — t, which implies that all terms must be 
even in a operators. 

• Point group operations (In one dimension limited to reflections around a point). 

Summing out the surfactant directions on site % gives rise to terms proportional to (1 — of) 
and involving the cr's on the neighboring sites. In one dimension, four terms of this type are 
consistent with the above symmetries: 

*o(l " o]) 



*l(l 



X 2 (l - of)0-j_i0-j + i 
^3(1 - of) of-i of+i 



(4) 



X 2 corresponds to L in the original Schick hamiltonian. Considering triplets of sites, the 
following four equations emerge from summing out the directional degree of freedom at the 
central cite, 

exp(-/5Kfr) = E exp(-/?W) , (5) 

T 

with (3 = l/k B T, T temperature, and k B the Boltzmann constant: 

©-©-©: X 2 + X 3 + 2X 1 +X Q = -k B T\n2 

- - © : -X 2 + X 3 + 2X X + X = -k B T ln(e 2/M + e" 2/M ) 

(6) 

©-©-©: X x +X = -k B T hi{e 0A + e' 0A ) 

- - : X = -k B Tln2 

Here the symbols 0, 0, and © represent a = — 1, a = 0, and a = +1, respectively. Solving 
the linear equations (^) is trivial: 

X = —k B T In 2 

X 1 = -A- k B T ln{l + e~ 2(3A ) (7) 
X 2 = A + h B T\n(l + e-^ A ) 

1 1 + e -W 

X 3 = A - -k B Tln(- 



2 v (l + e" 2 ^) 4y 

Note that the term proportional to X\ in ([|) splits into a single-site term and a pair term, 
and the term proportional to X splits into a constant and a single-site term. Thus, in one 
dimension the effective spin-1 hamiltonian is: 

H eS = E {-Ha, - (/i + X - 2X 1 )<t? - Ja i( x J+1 - (X + 2X 1 )<jffi +1 

(8) 

+ (1 - <7 2 )(X^_ 1( 7 i+1 + Xsal^A + X ) 
This hamiltonian has an extra three-spin interaction, X 3 crf_ 1 (l — of )of +1 , compared with 
the original Schick hamiltonian, and in higher dimensions even more interactions must be 
included. 

In the above expressions for the effective couplings (equations [7] and [J7p, the first term is 



the zero-temperature value, and the temperature dependence is isolated in the second term. 



With exception of the entropic single-site X coupling, all the effective couplings have a 
simple monotonous temperature dependence, and tend to zero with increasing temperature. 
Since the temperature behavior of the effective couplings in the CHS — ► Schick mapping 
is rather trivial, many qualitative features of the CHS model are probably preserved when 
one considers a spin-1 model with temperature independent couplings. Note, however, 
the nontrivial behavior of the nearest neighbor of a 2 coupling, (K + 2Xi), in the spin-1 
hamiltonian (equation H). For A > K/2 the coupling is negative at low temperatures, and 
changes sign with increasing temperature. 

The CHS model may be applied to the two component oil-surfactant (or water-surfact- 
ant) system (H — > ±00), and contains all the necessary ingredients to get a microemulsion 
like phase where water rich regions are separated by surfactant bilayers 0. In a study of 
this system using the Schick model [|TTJ, a four-site interaction was added to get the bilayer 
structure. The present calculation in one dimension indicates that it may be possible to 
model the bilayer structure with a spin-1 hamiltonian using two- and three-site interactions 
only. The important ingredient is probably the nontrivial of cr 2 coupling mentioned above. 



IV. MF THEORY OF THE EFFECTIVE HAMILTONIAN 

In this section standard mean-field theory is applied to the two-dimensional version of the 
effective hamiltonian of section |H|. In earlier work |3|,|4|||[^] mean-field theory was applied 
to the CHS hamiltonian ([I]) directly, and I expect that the present approach will improve 
on these results. 

The starting point of mean-field theory is to express the energy and entropy per spin 
as functionals of the magnetization, m = (a) ((. . .) represents thermal average), and sur- 
factant density, p s = (1 — cr 2 ), by assuming that the probabilities of the different states are 
uncorrelated at different sites. For homogeneous phases (all sites equivalent) the mean-field 
energy per site, e, is 

e = -Em - 2 Jm 2 - (1 - p s ) // - 2K(1 - p s f + 
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(x + 4X X (1 - p s ) + 2X 2 m 2 + 6X3(1 - Ps f + 4X 6 m 2 (1 - p s ) + 
4X 8 (1 - Ps f + X 9 m 4 + 2X 10 m 2 (l - p s f + X 12 (l - p s f) p s , (9) 



where the effective couplings, X n , are given by equation |37|. The entropy per site, s, takes 
the simple form 

/, 1-m-ps 1 -m-p s 1+m-ps l+m-p s 
s/k B = ln( ) ln( )-p s \n(p s ). (10) 

Note that when mean-field theory is applied to the CHS hamiltonian directly, the equiv- 
alents of equations |9| and are independent of the surfactant strength A. Hence the part 
of the phase diagram that pertains to homogeneous phases, ie disordered to oil/water co- 
existence transitions and tricritical point, is totally insensitive to the surfactant strength in 
that approximation. The present approach removes this artifact. 

The free energy at a given field, H, and chemical potential, //, f(/3,p:,H), is found by 
minimizing f(/3,p,,H;p s ,m) = e — (3~ l {s/k-Q) with respect to p s and m. Thus, the partial 
derivatives of / with respect to m and p s must vanish: 



dl 
dm 



-H - AJm + (4X 2 m + 4X 9 m 3 + 8X 6 m (1 - p s ) + 4X 10 m(l - p s ) 2 ) p s 



1 1-m-p 1 _l+m-p. 



dp s 



2p y 2 ' 2p y 2 
(11) 

pi + X + 4Xx (1 - 2p s ) + 2X 2 m 2 + 6X3 (1 - 3p s ) (1 - p s ) 

+ 4X 6 m 2 (1 - 2p s ) + 4X 8 (1 - Ps f (1 - Ap s ) + X 9 m 4 

+ 2X 10 m 2 (1 - 3p s ) (1 - p s ) + X 12 (l - p s f (1 - 5 Ps ) + 4X (1 - p s ) 

~ 2i5 ln( 2 } - W H 2 } + Hps) 

= (12) 

The m and p s that satisfy these equations will be the equilibrium values for m and p s . 
There may be several solutions, and the solution that gives the lowest free energy yields the 
equilibrium values. 
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By symmetry the critical line, separating the high-temperature disordered phase from 
the low-temperature oil- or water-rich phase, is found at H — m — 0, where equation |11| is 
trivially satisfied. Thus, only equation ( |T2| ) is important for the localization of the critical 
line. At the critical line the minimum changes into a maximum. The minimal eigenvalue of 
the matrix of second partial derivatives of / must therefore vanish. In the present case the 
matrix has a very simple form: the nondiagonal elements vanish, and the diagonal element 
(d 2 f /dpi) is always larger than the element (d 2 f /dm 2 ). Hence, the critical line is defined 
by the vanishing of the latter: 



-47 + \ + (4X 2 + 8X 6 (1 - Pa ) + 4X 10 (1 - p s f) Ps = 



(13) 



P(l-P.) 

This critical condition is independent of the chemical potential, p, and it is therefore, for 
given interaction strengths and temperature, an equation for p s . The equation is a fourth 
order polynomal equation, and of the four solutions, only one is physically acceptable (p s 
real, and < p s < 1). This solution is the critical density, p^ nt (T). The surfactant density 
and chemical potential are connected by equation [12] giving the critical chemical potential 



/i crit (T; p s ) = —AK (1 - p s ) -X - 4Yx (1 - 2p s ) - 6X 3 (1 - 3p s ) (1 - p s 

-4X 8 (1 - p s ) 2 (1 - Ap s ) - X 12 (l - Ps ) 3 (1 - 5p s ) 
, 1 i ( l ~ P g \ 



(14) 



The phase transition is not always continuous. For low temperatures it changes to first 
order at a tri critical point, below which the free energy may be lowered near the "critical" 
line, as defined by equations (0) and fli~4j) , by an infinitesimal displacement (m,5p s ) = 
(e, ae 2 ), ie a quadratic relation between 5p s and m (5p s ~ m 2 ). All odd partial derivatives 
with respect to m vanish by symmetry, and close to the critical line the free energy functional 
is then: 

1 ( d 3 f \ 2r , 1 (<)-.(' 
\dm 2 df 

I f0 2 f\ o 1 { d 3 f \ 1 ( d A f 

a H 



f { m, 5 Ps ) = /(0, 0) + \ (g) 5p 2 s + t j m>6p. + ^ 



/(0,0) 



2\dp 2 J a 2\dm 2 dp s 



24 \ dm 4 



J m 4 + 0{m b ) 
m 4 + C(m 6 ) (15) 
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The tricritical point is the point where, coming from higher temperatures, it is first possible 
to find an a such that the term in square brackets ([ ]) in this equation ( |T5D vanishes. This 
defines the equation 



( &f V ifd*f\ f&£ 

\dm 2 dp s J 3 \dm 4 J \dp 2 s / 



(16) 



with 



d 3 f 



dm?dp s ) (3{l-p s f 

[d^J ~ (3(1- Ps ) 
fd 2 f\ 1 



+ 4X 2 + 8X 6 (1 - 2p s ) + 4X 10 (1 - 3p s ) (1 - p s ) (17) 
3 + 24X 9 p s (18) 
- AK - 8X1 - 12X 3 (2 - 3p s ) 



dpi/ (3(l-p s )ps 

-24X 8 (1 - p.) (1 - 2 Ps ) - 4X 12 (1 - Ps f (2 - 5p s ) (19) 



to be satisfied at the tricritical point in addition to equation [13]. The tricritical density, 
p S 3, and temperature, T 3 , as a function of surfactant strength, A, and of the parameter K 
are given in figure ^. The Lifshitz point is calculated in section |V|. Note that in earlier 
mean-field calculations j|j7j the position of the tricritical point was independent of A, and 
the phase diagram was always of type C for large A. In the present calculation this type 
of phase diagram is found in an intermediate A range only, and only for K w 0. For a 
sufficiently strong surfactant the phase diagram is always of type B. 



V. CORRELATIONS IN HIGH-TEMPERATURE PHASE FROM MF-THEORY 

The connection between the free energy functional and the spatial correlations is textbook 
material ||22|| , and has previously been applied to mean- field theories of microemulsion models 

ISMS 

While the equilibrium values of the fields (m and p s ) are the same at all sites, the 
fluctuations are not homogeneous. Expanding to lowest (second) order in the fluctuations, 
5m and Sp s , the free energy takes the form 
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F = F (m,p s ) + J2<PnHd)(f : 



n+d 



n,d 



(20) 



with 



$(d) 



^$ mm (d) $ mps (d) N 
v $ Psm (d) $ PsPs (d) y 



(21) 



and 



n = (Sm n ,5p sn ) 

The Fourier transformed correlation function can then be expressed as 

(27T) 2 



(22) 



H(k) = 



-$- x (k) 



(23) 



In the present short-range model $(k) is a finite sum over near neighbor terms. We are 
primarily interested in the zero field model (H = 0) where m — 0. Here $(k) is diagonal 
(non-diagonal elements are proportional to m) with the elements: 

1 



—2 J (cos(/c x ) + cos(fcy)) 

+ (2X 2 + 4X 6 (1 - p s ) + 2X 10 (1 - Ps ) 2 ) Ps (cos(2k x ) + cos(2k y )) 



(24) 



(25) 



P(l-Ps)p 8 

- ( 2X + 4Xi + 12X 3 (1 - p s ) + 

12X 8 (1 - p s f + 4X 12 (1 - p s f ) (cos(A; x ) + cos(Ag) 
+ (2X3 + 4X 8 (1 - p s ) + 2X 12 (1 - p s f) p s (cos(2k x ) + cos(2Ag) 
+ (4X 3 + 8X12 (1 - p s )) p s {cos(k x - ky) + cos(k x + k y )) 

A detailed discussion of the form of these correlations, and comparison with experimental 
and Monte Carlo scattering functions, and other theoretical work, is deferred to a separate 
paper. 



12 



A. The Lifshitz line 



The competition between the cos(/c)-term (proportional to the water-water interaction, 
J) and the cos ( 2 k) -term (proportional to the surfactant density, p s ) in $ mm (k) (equation 
|2"4"|) is the origin of a maximum in the correlation function, H mm (k), at k 7^ in parts of 
the phase diagram. This maximum is an experimental characteristic of microemulsions, and 
its location determines the characteristic length scale for the microemulsion structure. The 
maximum is located on the diagonal (k x = k y = k maiX ): 

cos(/c max ) = t - — (26) 

(4X 2 + 8X 6 (1 - p.) + 4X 10 (1 - psf) p s 

The Lifshitz line separates the region of the phase diagram where the correlation function 
has a maximum at k / (microemulsion) from the region with a maximum at k = 
(ordinary fluid). Note that no phase transition occurs at this line. The Lifshitz line satisfies 
the equation 

cos(fc max ) = 1 . (27) 



B. Phase transition from disordered to incommensurate 



Critical behavior is recognized by divergencies in the correlations, i.e. by zeros in $(k). 



The disorder-oil/water critical line (equation [13]) is recovered as $ mm (0) = 0, while di- 
vergencies at k 7^ signal a possible continuous transition into an incommensurate phase 
(ordered structured phase where the order is incommensurate with the underlying lattice). 
In this section I assume that the transition is continuous. There is no real justification for 
this assumption, and the transition will most probably be first order in at least part of the 
phase diagram 

The transition will take place at k = (fc max , ^max) (equation |26]), and for a diagonal k the 
equation $ mm (k) = is quadratic in cos(fc). In different regions of the phase diagram the 
solutions to this equation fall into one of three classes: 
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1) Two complex conjugate solutions: No maximum at k 7^ 0. 

2) Two real solutions: $ mm (k) < for some k. This is unphysical and the system 
must be in an ordered phase. 

3) A single real solution: Divergence at k 7^ 0, i.e. critical line. 

Using the criterion for class 3 above, the line of continuous phase transitions into an 
incommensurate phase satisfies the equation 

b 2 - 8c(a - c) = (28) 

with 



P(l -Ps) 
b = 4J 

J 



Ap s (X 2 + 2X 6 (1 - Ps ) + X 10 (l - Ps 



cos(fc 

1] 



VI. MONTE CARLO SIMULATION 

The model has been simulated using a standard Metropolis Monte Carlo method. The 
lattice is divided into three sublattices to enable vectorization, in a manner very similar to 
the partition used by Wansleben et al ||2^j2^j , but no multispin coding is used. The random 



number generator employed is a shift-register generator with a very long period (2 607 — 1), 
which is never exhausted in a MC run (For a discussion of random generators of this type, see 
the paper of Compagner and Hoogland |26] and references therein). The program generates 



~ 10 6 spin-flip trials per second on a CRAY X-MP. To simplify (and speed up) the program, 
I have only simulated the K = 0, version of the hamiltonian ([J). 

Detailed Monte Carlo studies have been performed for the K = 0, A = 3J model only. 

Standard finite size scaling techniques was used to locate the critical line pTfl . This 
includes scaling of the magnetization (m), susceptibility (X m ), and fourth order cumulant, 

{U m ): 
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m 



(29) 



X m = (rn 2 ) - (my (30) 
U m = 1 - (31) 

L is the linear size of the lattice. Since all the critical exponents are exactly known (the 
order-disorder transition is in the universality class of the two-dimensional Ising model) 
there is only one parameter, the critical temperature, T c , to be fitted in all these scaling 
laws: 

m(r) = L 1/8 m(r) 
Xjr) = L- 7 ^X 



m \ I 



U m (r) = U m {r) (32) 

'T-T n 



T 



T, 



c 



L 



As is evident from figure ^ finite size scaling works well for all scaling functions, with a 
unique critical temperature. 

The next question to be answered is whether the phase transition changes into first order 
at a tricritical point, as is the case in the mean-field phase diagrams calculated above. The 
Monte Carlo results show pronounced hysteresis at low temperatures, but this is definitely 
not conclusive, and especially not so in this particular model, since the tricritical temperature 
is very low. More convincing is the fact that the susceptibility, \ , of the noncritical 
surfactant density, p s 

Xs = (PsPs) - (Psf (34) 

at the critical line obeys finite size scaling using the tricritical exponents relevant for the 
model PS] (Universality class of a one component order parameter in two dimensions, ie 
tricritical spin-1 Ising model): 
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X s (r) = L-^Xs 



or 

T — T a 



(35) 



T 3 

with z/ t = 5/9 and j t = 37/36. Tricritical finite size plots are shown in figure [5]. The resulting 
tricritical point for K = and A = 3J is k B T 3 /J = 0.52 ± 0.03 and /i 3 /J = 3.15 ± 0.1. 

The Monte Carlo estimates for points along the line of continuous oil/water — > disordered 
transitions, and of the tricritical point, for A = 3.0J and K = are included in figure || 
The present calculation is inconsistent with earlier Monte Carlo calculations by Laradji et 
al ||, which gave critical temperature estimates that are about a factor 1.4 higher than the 
present estimates. I have no explanation for the discrepancy, but tend to believe that the 
present estimates are better than the estimates of Laradji et al. The reason for this is that in 
the present calculation the critical temperature is always lower than the pure Ising (fi = oo) 
value. This is not the case in the calculations of Laradji et al. Normally the presence of 



surfactants will counteract ferromagnetic order, as it does in the A = model , and it is 
difficult to see how it could promote such order. 

It was not possible to get conclusive results on the nature of the d — > i (disordered- to 
ordered structured phase) transition. I will therefore not rule out completely the possibility 
that a first-order d — > i line crosses the critical line at a higher temperature than the tricritical 
point found above. The model becomes very difficult to study with the present algorithm 
close to the tricritical point. In the order of > 10 6 whole lattice sweeps were necessary 
to get reasonable statistics at the points simulated, and the finite size plots in figure |] 
represents ~ 20 hours of Cray CPU-time. The dynamics grows slower as the tricritical point 
is approached. The reason for the slow dynamics seems to be that the activity is limited 
to oil /surfactant /water interfaces while most of the system lies within homogeneous oil or 
water regions. 
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VII. CONCLUSION 



Representative phase diagrams are shown in figures |6] and 0. Except for weak surfactant, 
the phase diagrams are all of type B, and include an oil/ water /mi croemulsion coexistence 



region. This typical behavior was also found in studies of the Schick model [10,18,19 



and this supports the claim in section [TTT| that this simpler model preserves the important 
physical properties of the CHS model. The previous simple mean-field result [|,§],[7| that 
the phase diagram is always of type C for strong surfactant, seems to be an artifact of the 
approximation employed. 

Results from Monte Carlo simulations are consistent with the phase diagram topologies 
found from mean-field theory, but the results are not conclusive. The quantitative error in 
the location of the mean-field critical and tricritical temperatures is of the order expected in 
two dimensions. Note that the tricritical temperature found here in the simulations on the 
K = 0, A = 3 J model, k-^T^/J = 0.52 ± 0.03 is lower than the tricritical temperature for 



the K — 0, A = model |30 |, k-oT^j J = 0.610 ± 0.005. This confirms the mean-field result 
that T3 initially decreases with larger A (see figure |3|). The critical and tricritical surfactant 
densities found in the Monte Carlo simulations are much lower than the corresponding mean- 
field densities. The amount of surfactant necessary to make oil and water miscible is thus 
much lower than what mean-field theory predicts. For the K = 0, A = 3J model the 
tricritical density is Pg 3 ~ 0.06 compared with the mean- field result pf£ = 0.2, and the 
A independent K = tricritical density p^ 3 F ' = 2/3 of mean- field theory applied to the 
CHS hamiltonian directly is totally off the mark. The Monte Carlo results are compatible 
with experimental microemulsion systems where the amount of surfactant is typically a few 
percent. Notice also (figure |6|) that the MC critical line is almost vertical (constant p s ) over 
a wide range of temperatures near the tricritical point. This feature is not shared by the 
mean-field critical line which has an almost constant slope for all temperatures. 

The present study strengthen previous results on the same model that the decisive mi- 
croscopic ingredient for microemulsion behavior is the amphiphilic nature of the surfactant. 
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In two dimensions the new effective interactions resulting from the summing out of the 
surfactant directions at a site i involve the factor (1 — of) and all combinations of a and a 2 
on i's next neighbor sites that are compatible with the symmetries of the CHS hamiltonian. 
Labeling the central site o and numbering its nearest neighbor spins o\ . . . 04 in the clockwise 
direction, these combinations are: 

*o(l 

X x {\ - o 2 )(o 2 + aj + al + a\) 

X 2 (l- Oq)(oio 3 + o 2 o 4 ) 

X 3 (l - o 2 )(o 1 2 o 3 2 + o 2 2 o 4 2 ) 

X 4 (l - o 2 )(oio 2 + 02O3 + 03CT4 + 0401) 

X 5 (l - cr 2 )(a 2 a 2 2 + o 2 o 2 + a\a\ + a\a\) 

X 6 (l - a^ia^a 2 + o 4 2 ) + a 2 a A (a\ + of)) (36) 
X 7 (l - o- 2 )(oio- 2 (o 2 + a\) + a 2 a 3 (a 2 + o 2 ) + a 3 a 4 (a 2 + a 2 ) + o±o x {p\ + o 2 )) 
X 8 (l - a 2 )(a 2 a 2 3 (a 2 2 + o 2 ) + o 2 o 2 (o 2 + o 2 )) 

-^9(1 — 0"o)o"lO"20"3°"4 

-^lo(l - O'o)! "! "^^ "^ + 01020304) 

, 2 \ / 222 222 2 2 \ 

Xn(l — O o )(01020304 + 01020,304 + 01020304 + 010 2 3 04) 

X 12 (l - (y 2 ^)a\a 2 2 a\a\ 
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Considering configurations of four spins around a central surfactant, the summing out of 
the surfactant directions leads to a set of linear equations for the effective couplings, X n . 
Solving these equations gives: 

X = -k B T\n4 

X l = -A + 2k B TH T ^ x ) 

X 2 = A-A; B Tln( rT ^- w ) 

8 (l + e -2/M\ 
X 3 = A - k B T\n( — /) 

X 4 = 



X5 — x 3 



11 4(1 

•Vi; = —-A + -k B T\n(— — j-j— — —777-7 



2 2 " (l + e -2/3Ay 

X 7 = (37) 



, 6 



3 1 1024(l + e~ 2/M ) 

X8 = ~2 A + 2 kBTH (l + e -^-e-^)(l + e~^r' 

1 1 8 (l + e- 4 ^) 
X g — -A — -k B T ln(-^ 1) 

9 2 4 v (1 + e" 2 ^) 4 



4 



1 1 32 (l + e -2/3A_ e -/3A\ 

Xi = -A - —k B T ln( ^ f T ) 

x n = o 

5 1 2 35 (l + e~^ A ) (l + e~ 2 ^ 2t ' 



X 12 = -A - —k B T ln( i ^ ^-gj) 
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FIGURES 

FIG. 1. Types of phase diagrams that is expected in the CHS-model [§,0]. T is temperature, and 
[i is the surfactant chemical potential, type A, (a), is relevant for a weak surfactant: two regions, 
disordered (d) and oil/water coexistence (o-w) is separated by a continuous phase transition (full 
line) which changes into first order (dashed line) at a tricritical point. With a strong surfactant, 
the phase diagram is either of type B, (b), or type C, (c). Both these diagrams show four regions: 
disordered (d), microemulsion (m), oil/water coexistence (o-w), and ordered structured phases (i), 
which may be incommensurate with the underlying lattice. 

FIG. 2. Ground states, and T = phase diagram of CHS-model in two dimensions. The model 
has three different ground states: All oil or water (F), all surfactant (S), layered (L), and tubular 
(T). The layered and tubular states, correspond to ordered structured phases (i) in the T / 
phase diagrams of figures |l[ || and ^. 

FIG. 3. Tricritical density, p s s, (top) and temperature, T3, (bottom) as a function of surfactant 
strength, A, with K = 0.0. . . 1.0 (Full lines). The dashed lines give the position of the Lifshitz 
points in regions of parameter space where the phase diagram is of type C (fig. |l|). 

FIG. 4. Finite size scaling plots for the magnetization, fh = L^m, fourth order cumu- 
lant, U m = U m , and susceptibility, ^ m = L -7 / 4 ^ m as functions of scaled temperature, 
r = L(T - T c )/T c . At K = 0, A = 3.0 J and p = 5.0 J, with T c = 1.265 J. 

FIG. 5. Tricritical scaling: Finite size scaling plots for the surfactant density susceptibility, ^ , 
at criticality as a function of scaled chemical potential (top) and scaled temperature (bottom). At 
K = and A = 3.0J with ^3 = 3.15 and T 3 = 0.52. 
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FIG. 6. Phase diagrams for K = and three different surfactant strengths. The lines are the 
results of mean-field theory on the effective spin-1 hamiltonian. Full lines are continuous- and 
dashed lines first-order transitions. The dotted line is the Lifshitz line. Note that the line for the 
d — > i transition is found under the assumption that the phase transition is continuous, that the 
region denoted by i (incommensurate) is probably divided into several different phases separated by 
first-order phase transitions, and that this region is separated from the o-w region by a first-order 
transition This is not indicated in the figure. For comparison the (A independent) localization 
of the tricritical point from the application of mean-field theory on the CHS-hamiltonian directly is 
indicated (with a +). Squares are Monte Carlo results for the d — > o-w transition (A = 3.0 J only). 
Error bars for the critical temperature are not shown, but are much smaller than the symbols, while 
the error bars for the critical surfactant density are shown. The MC tricritical point is indicated 
by a cross surrounded by an error rectangle (A = 3.0 J and top row only). 

FIG. 7. Phase diagrams for K = J and three different surfactant strengths. The notation 
follows figure |(| 
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